/* 
 * File:   integracion.h
 * Author: jonk
 *
 * Created on 15 de enero de 2012, 23:43
 */

#ifndef INTEGRACION_H
#define	INTEGRACION_H

double integracion(particula *part, parametros param)
{
    double epoten;
    
    epoten= interaccion(part,  param);
    
    for(int i=0;i<param.npart;i++)
    {
        for(int j=1;j<3;j++)
        {
            part[i].pos[j]+=part[i].force[j]/param.zeta*param.dt;
        }
        part[i].pos[0]+=(part[i].force[0]/param.zeta)*param.dt;//+part[i].pos[2]*param.gamma0*param.omega*cos(param.omega*param.t))*param.dt;
    }
    
    return epoten;
}

#endif	/* INTEGRACION_H */

